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Abstract 

The multiple particle final states in modern data from 4ir detector ex- 
periments offer a wealth of information about higher hadronic resonances 
and hadron interactions. However, it requires careful analysis to extract 
model independent results from those data. 

1 Introduction 

The handling of i-exchanges in scattering experiments requires serious consider- 
ation. It is the first and dominant process involving multi-particle states, which 
is relevant for the study of hadronic resonances at energies of the second reso- 
nance region (1.5 GeV and above). In most approaches three and more particle 
states are handled as effective two-particle states. The state with one nucleon 
and two pions is often modeled as a rho-nucleon or a sigma-nucleon state. How- 
ever, if the higher resonances are to be understood in a model-independent way, 
the three-particle states, and the non-resonant production processes, such as 
the i-exchange processes should be handled in their full glory. Fundamentally, 
this is not a big problem; the theory is well understood. However, practically it 
is a different story. Central to the computational problems lie the singularities 
associated with the three-particle states going on-shell. It should be noted that 
in a simplistic Yukawa treatment of the i-exchange, ignoring the energy depen- 
dence of the exchange diagram, such a problem does not occur. However, if the 
three-particle states and the singularities are treated seriously, they should be 
regulated and integrated over which requires three scales: the regulation scale 
e <C the integration scale AE the physical scale /i. Every problem is big in 
numerical sense, because the handling of triple scales. Furthermore, if one works 
with momentum variables, the singularities form curves in the kincmatical do- 
main, so extra care is required to locate the singularities and extrema on the 
curve with respect to the chosen grid. In terms of energy variables the position 
of the singularity is obvious; when the scattering energy equals the three-particle 
energy, however, the numerical integration still requires great care. In recent 
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years some methods, such as the Lorentz transform, have been developed, which 
essentially takes one away from the real axis and the singularities, however, they 
are rather complicated. 

In this paper another approach is investigated. Central is the expansion 
of the amplitudes and relevant function in an orthogonal basis with known 
dispersion integrals, which absolves one from doing any singular integration. In 
terms of bookkeeping it does require a serious effort, as one should keep track 
of the real and imaginary parts of both the two-particle and the three-particle 
thresholds and amplitudes. The threshold of the three-particle state depends 
on the two-particle states from which it originates; if the two-particle state has 
zero momentum, i.e., E = m a + mb, the three-particle state it produces can have 
zero momentum for all particles too, and its threshold lies at E = m a + nib + fi, 
while if the two-particle state has a higher momentum, the threshold of the 
three-particle states it can produces will lie appropriately higher. Furthermore, 
given the energies of the initial and final two-particle state and w/, there is 
an upper bound to the energy of the three-particle state it can produce in the 
i-exchange where the exchange particle is emitted from particle o: 

^max/min = y/fc± kf) 2 + /I 2 + yj kf + m 2 + ^ k 2 + m\ , (1) 

where k is the momentum associated with the energy u: 
= y/(a> 2 - m 2 a - m 2 ) 2 - 4m 2 ^f 

Eventually, the kinematical restrictions on a two-particle state in a given 
partial wave will lead to a set of three-particle states, which can be labeled by 
the angle between the momentum of the exchange particle and the momentum 
of the two-particle state. This integration can be performed analytically and 
leads in the case of a two-particle s-wave to a second-order transition amplitude: 

T(°)(£;,a; < ,a; / ) = ^-Iog 

where g is the coupling constant, and numerical factors are ignored. The three- 
particle state in now implicitly defined through the imaginary part of the two- 
particle to two-particle transition, which saves one from constructing a basis for 
the three-particle states. As every three-particle state follows from the emission 
of an exchange particle from the two-particle state, the three-particle amplitude 
is given by the two-particle amplitude times the imaginary part of the elemen- 
tary transition amplitude T^°\ The two-particle amplitude is determined from 
solving the Lipmann-Schwingcr equation with the elementary transition ampli- 
tudes like in the kernel. 

The real part of the amplitude is dominant. The leading order imaginary 
part is the first order process given by the imaginary part of However, as 
with most problems, the qualitative changes in the results through the imaginary 
parts, i.e., on-shell three-particle states, can only be understood properly if they 
are included in the calculation. 
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Making a separable expansion of the amplitude lies at the basis of an effec- 
tive method to sum the processes to all order to yield a unitary and analytic 
transition amplitude: 

T(°\E,u u u Jf )=^a lm 4 E \u Ji )4>£\u; f ) , (4) 

Ira 

( E) 

where one should take care the the basis (jr a (to) is orthogonal (or band-diagonal, 
such as splines) for a well defined expansion, and it has known dispersion inte- 
grals: 

The real part of the amplitude spans the full kinematical domain m a + rrib < uj, 
hence an expansion in £ = (to 2 — (m a + mi,) 2 )/^ 2 e [0, 1] can serve as a basis 

^) = V4TT2 (ma +f b) P ; (2£-l) , (6) 

for products of which a closed form dispersion integral can be calculated. There- 
fore, for the expansion of the real part of the suffix is not necessary. 

However, the imaginary part of is only non-zero in a restricted but 
infinite domain of E ® u>i ® u>f, when the scattering energy E lies in between 
S m i n (wj,w/) and E maK (ui,ujf). The domain in the Wj ® uij space given a scat- 
tering energy E is the energy equivalent of the smooth triangular kinematical 
region of a Dalitz plot, where the each of the corners is sharper if the energy 
is larger compared to the mass m a , nib, or \i. As the energies E, Wj, and u)j 
increase, the domain approaches the triangular area: 

E < uii + ujf and E > Ui and E > ujj . (7) 

This correlation between the variables does not allow for a separable expansion 
of the imaginary part of . 

It seems that a useful parametrization of Dalitz plot amplitudes and it gen- 
eralization to more particles and more dimensions does not yet exist. In the 
massless case, the domain is a simplex on which Appell polynomials and gen- 
eralized mappings of orthogonal polynomials from the n-ball to the n-simplex 
form a multi-dimensional basis. In the massive case a system of a weight func- 
tion that serves also as domain definition times polynomials Pi in the single 
particle energies e of the three or more particle state might serve many needs: 

</>! = (E - E cdgc ) (E- E cdgc f ff> (£!,...,€„_!) , (8) 

such that 

dei ■ ■ ■ de n -i<j)i<j)j = 5 U . (9) 
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For each case of different masses (ttttN, irpN, etc) the coefficients of the orthog- 
onal polynomials Pj must be determined. Furthermore, the dispersion integrals 
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Figure 1 : The domains where the imaginary part ofT(°) is nonzero for increasing 
scattering energies (m a = nib = 10/u). 

must be calculated for both: 

— / du> f^^ 3 and — / duj-J^- — . (10) 
ir J E - lo ir J E -u> 

However, after this preliminary work, the calculation of analytic and unitarity 
transition amplitudes that incorporate resonant and i-exchange contributions, 
and three-particle final states will be a straighforward problem. This work is 
under investigation. 
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